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Abstract 

A new model describing the dark sector of the universe is established. The model 
involves Bose-Einstein condensate (BEC) as dark energy (DE) and an excited state above 
it as dark matter (DM). The condensate is assumed to have a negative pressure and is 
embodied as an exotic fluid with Chaplygin equation of state. Excitations are described 
as a quasiparticle gas. It is shown that the model is not in disagreement with the current 
observations of the cosmic acceleration. The model predicts increase of the effective 
cosmological constant and a complete disappearance of the matter at the far future. 
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1 Introduction 



The energy content of the Universe is a fundamental issue in cosmology. Observational data 
are evidence of accelerating flat Friedmann- Robertson- Walker model, constituted of about 1/4 
of baryonic and dark matter and about 3/4 of a dark energy component. The DM content was 
originally inferred from spiral galactic rotation curves and then was supported by gravitational 
lensing and cosmic microwave background observations. 

The essential feature of DE is that its pressure must be negative to reproduce the present 
accelerated cosmic expansion. The simplest DE model, the cosmological constant, is indeed 
the vacuum energy with the equation of state p = —p. The models for which p < —p has been 
denoted phantom energy, and possesses peculiar properties, such as, an infinitely increasing 
energy density [1], negative temperatures [2], and the violation of the null energy condition. A 
number of models, such as quintessence [3] and k-essence [4], are based on scalar field theories. 
These models are parameterized by an equation of state p < — p/3. For a recent review of DE 
models and references see [5]. 

An alternative model is that of the Chaplygin gas, also denoted as quartessence, based on a 
negative pressure fluid, which is inversely proportional to the energy density [6]. The equation 
of state representing the generahzed Chaplygin gas (GCG) is given by pch = ~^/Pch '^it^ 
positive constants A and a (0 < ck < 1) [7]. An attractive feature of these models, is that at 
early times, the energy density behaves as a matter, pch oc a~^, where a is the scale factor, 
and as a cosmological constant at a later stage, pch =const. It is also suggested that at an 
intermediate stage the energy density pch consists of both vacuum and soft matter (matter 
with the equation of state p = ap) contributions. This is favorable to use the GCG model for 
a DE and DM unification [7, 8, 9]. 

Some different approach to the same problem is realized in this work. The feature of the 
scalar fleld used in the model [8] is spontaneous symmetry breaking that leads to a nonzero 
expectation value of the field. In other words, Bose-Einstein condensation of the scalar field 
into the state with zero momentum takes place. Similar effect holds in the theories of super- 
conductivity and superfluidity. At zero temperature superfluid is in its ground state. If T 7^ 
particle-like excitations arises above the ground state. In this case the system can be divided 
into a background superfluid condensate and a quasiparticle gas. This separation naturally 
leads to two-fluid dynamics in which the condensate is said to be a superfluid component and 
the quasiparticle gas forms a normal component. 

The model proposed in this letter represents the dark sector of the universe as a superfluid 
where the superfluid condensate is considered as DE and the normal component is interpreted 
as DM. To provide the accelerated expansion the potential of the scalar field must have specific 
form and entail a negative pressure of the superfluid background. 

We base our analysis on the action 
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where Lagrangian £ associated with a generahzed hydrodynamic pressure function depends 
only on one variable if we consider pure condensate, and on three variables when we include 
the excitation gas. 



2 Condensate in WKB-approximation 

Consider a system with spontaneously breaking of C/(l) symmetry described by a complex scalar 
field with nonzero Gibbs expectation 0(x) = (■0), the condensate wave function. Quantum 
fluctuations -0' = -0 — can be considered in terms of quasiparticles. 

If the quasiparticle gas is dilute or interaction between the quasiparticles and the condensate 
is weak then influence of the elementary excitations on the ground state can be neglected. In 
this case the condensate is described by the Lagrangian 

£ = a^0*a>-v^(0». (2) 

It is useful for further study to represent the Lagrangian (2) in an equivalent hydrodynamic 
form. For this purpose we write the condensate wave function in terms of modulus and phase: 

^=^.-« (3) 
Substituting (3) into the field equation, the real and imaginary parts yield respectively 

V.VV + a (^2^ - d.xd''^ = 0, (4) 

and 

V.(a29^x) = 0. (5) 
Equation (5) is a conservation law for the 4-current 

3, = i {4>*dA - d,4>*4>) = a'^d^x- (6) 

The gradient of the condensate phase is a superfluid momentum which can be written in terms 
of a unit 4- vector V^: 

duX = fJ'u = /iK, (7) 
where is a chemical potential. In the present context ji, is the particle current, therefore 

where ric is a particle density of BEC. 

If the modulus a varies slower than the phase x, d^cr < crdi^x, then we can neglect derivatives 
di,a that corresponds to the WKB expansion of the condensate wave function up to the first 
order. 
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In this approximation equation (4) takes the form 



and the particle density is 



The energy-momentum tensor is found as variation of the Lagrangian with respect to the 
metric: 

where the pressure and the energy density 

depend on only one variable ^. 

For various contexts wc can apply a rather different kind of potential. In our case the 
condensate is assumed to possess a negative pressure. For this reason let us take the potential 

1^(0*0) = M f ^ + A y (13) 



^ A 

In accordance with equations (10) and (12) one finds that 

2Xii 2M 
n-c — — Pc = — (14) 

Vl-A/xVM yi-A//VM ^ ^ 

and the function of the generalized pressure has the form 

P{l^) = Pc = -2M^1 - A/iVM (15) 

For further purpose it is more rational to eliminate the chemical potential jj, and to express the 
hydrodynamic quantities via the energy density: 

AX / 4M^ 

^'^-\mVp'c-^M^ pe = , (16) 

V -IV-l Pc 

and the adiabatic speed of sound is 



2 clpc 4M 

= = (17) 



c 



The equation of state (16) is uniquely proper to Chaplygin gas suggested recently by Kamen- 
shchik et al. [6] as an alternative to quintessence and developed by a number of authors for 
description of the dark sector of the universe [7, 8] . 

In contrast to these works where pressure of Chaplygin gas is formed by both DE and DM, 
this model implies that the equation of state (16) concerns with only BEC which is interpreted 
as DE. 
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3 Relativistic superfluid dynamics 



An efficient approach to description of the excited state is two-fluid hydrodynamics. This theory 
does not depend on details of microscopic structure of the quantum liquid and exploits effective 
macroscopic quantities. In the theory there exist two independent flows, the coherent motion 
of the ground state named a superfluid component, and a normal component produced by the 
quasiparticle gas. For this reason it is necessary to increase the number of independent variables 
in the generalized pressure (15) from one to three [10, 11]. They correspond to three scalar 
invariants which can be constructed from the pair of independent vectors, namely superfluid /Xq, 
and thermal Oa momentum covectors so that the general variation of the generalized pressure 
in a fixed background is 

5P = n^'diia + s^'SOa. (18) 

The coefficients n" and s" arc to be interpreted as particle number and entropy currents corre- 
spondingly. By virtue of its invariance the pressure is given as a function of three independent 
variables, Ji = l/Xa/^", h = fJ^aO"', h = \OaO°'. Taking the derivatives of the pressure, one finds 

dP dP dP dP 

As soon as the generalized pressure is the Lagrangian density in the action (1) its variation 
with respect to the metric gives the energy-momentum tensor 

dP dP dP 

Tap = Ql^l^oclJ'p + -Qj-JylJ-aOfi + 9alJ.fi) + gj^^ocdp " -Pg'a/J- (20) 

Instead of the thermal momentum 6a let us introduce an inverse temperature vector /3" = 
s"/ (s^Op) which we use as the independent vector together with the superfluid momentum ^a 
since they are comoving to the excitation gas and the condensate respectively. Corresponding 
unit 4-velocities are 

C/° = /- , V = ^ (21) 

In place of the scalars Ii, I2, h we use new three invariants, a chemical potential /i — 
^J~ljiFiIp, scalar 7 = VaU°' associated with the relative motion of the components, and inverse 
temperature with respect to the reference frame comoving to the excitation gas /3 = \f^^^p. 

Using (19) and (21) the energy-momentum tensor and the particle number current are 
readily represented as 

n° = nc\/" + nn?/'' (22) 
Tap = IJn,VaVp + WnUaUp-Pgap, (23) 

Relations between the macroscopic quantities involved in equations (22) and (23) are of 
a quite general form. More detail information about them can be obtained from statistical 
description of the elementary excitations. The quasiparticle energy spectrum has a significant 
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nonlinear dispersion at high energy, and therefore completely relativistic description has been 
carried out only for a low energy excitations, phonons [12, 13]. Based on the relativistic kinetic 
theory of the phonon gas [13] we in particular can obtain 

^ = (1 - ciW., (24) 

when phonons prevail over another sorts of quasiparticles. 

Let us assume that the generalized pressure function is separated as follows: 

P(//, /3, 7) = p,{ix) + p„(//, /3, 7). (25) 

Equation (25) is best suited for the phonon gas since it describes neglect of excitation influence 
on the ground state. This ansatz retains equations (16) and (17) valid for the condensate in 
the framework of two-fluid dynamics. 

Let us also suppose that the pressure of the normal component depends on temperature 
by the power law, pn oc T"^+^. This assumption is a generalization of the dependence pn oc 
obtained by Carter and Langlois [12] for the equilibrium distribution of the phonon gas followed 
by the relation — c^W^/A. In our case it transforms to 

Equation (26) is a kind of barotropic equation of state and z/ is a properly polytropic index. 
Its value governs so-called second sound speed C2, the sound speed in the excitation gas. In the 
limit of low temperatures when the quasiparticle contribution becomes small C2 — )■ Cg/ yV, and 
in the case of the equilibrium phonon gas it coincides with the result obtained in [12]. 

We restrict our consideration to the equation of state (26) for the normal component situated 
between the dust one and the stiff one. It is evident from equation (26) that this constraint 
implies u >\. 



4 Universe with BEC 

4.1 Equations of motion 

The cosmic medium is now regarded as a matter which particularly is in the BEC state and 
its particle number current and energy-momentum tensor have the form (22) and (23). We 
assume that the self-dependent condensate ansatz (25) holds and the superfluid background 
obeys the equation of state (16) and the excited state is described by the relations (24) and 
(26). This means that we ignore nonlinear high-energy part of the quasiparticle spectrum. 
Although these restrictions render the model incomplete they essentially simplify the evolution 
equations retained simultaneously the key features of the model. 
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Let us consider a homogeneous and isotropic spatially flat universe. In this case the super- 
fluid and normal velocities are equal and thus 7 = 1. Einstein equations then reduce to 

■ 2 

3^ = STTGptot, -6- = 87rG(3ptot + Act), (27) 

where ptot consists of the condensate density Pc and the normal one pn = W^n — Pn that are 
interpretable as DE and DM densities respectively, and ptot = Pc + Pn- In accordance with the 
integrability conditions of Einstein equations we require local energy-momentum conservation 
V^T'^'^ = that yields 

Ptot + 3-(ptot + Ptot) = 0. (28) 
a 

The interaction between DE and DM is implicitly included in equation (28) and also in particle 
number conservation Vi^n^ = that leads to 

ritot + 3-ntot ^ =^ nc + nn= no = const. (29) 

This approach distinguishes the present model from [15] where the rate of the transition between 
ground and excited states F is explicitly used as an interaction factor and equation (28) breaks 
down into separated balance equations for DE and DM. In [16] the similar splitting is applied 
for interaction between Chaplygin gas (it is regarded as DE) and CDM. 

Taking into account the expressions (24)-(26) and (29) we reduce equations (27) and (28) 
to following two dimensionless equations: 



a 



2 



1 k ( fp y/ p^ — 1 



3«A+. *^V^fl-Af^-,34^|H0, (31) 



a 



where p — pc/2M, and k — no/2y/XM. The dimensionless time variable t' is connected with 
real time t as t' — y/ldnGMt. 

4.2 Exact solution 

In the formal limit z/ — > 00 equations (30) and (31) are solved analytically. As obvious from 
(26) the quasiparticle pressure is neglected and DM behaves as dust-like matter. In this case 
Eq. (31) yields 

and the scale factor varies in according to the integral 

t'^ f (33) 
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The parameter k gives an initial normalized total particle number density and h asso- 
ciates with an initial particle number density for the normal component. More precisely, 

— (no/nn(0) — If 6 = 0, the normal component is unavailable and the evolution 

follows the scenario proposed in [6] since the ground state obeys the equation of state (16). In 
the case of 6 7^ DM is governed by the law 



At the beginning stage (i.e. for small a) the total energy density is approximated by ptot c< 
that corresponds to a universe dominated by dust-like matter. The same behavior is 
a feature of Chaplygin gas [6] but even though in this model the condensate has the same 
equation of state, such dependence is due to the normal component. 

At the late stage (i.e. for large a) ptot 1- Separating now DE and DM contributions one 
finds the subleading terms are 

Pc ~ 1 + (35) 

Pn ~ ^, (36) 

whereas the scale factor time evolution corresponds to de Sitter spacetime, namely, a oc e*'/^. 
At this stage the fluid is almost in the ground state so that the condensate wave function is 
close to the potential minimum. The behavior is similar to GCG [7]: expressions (35) and 
(36) imply that the system evolves as a mixture of the cosmological constant and the dust-like 
matter. Note, the asymptotic formula (35) is valid for any value of u. 



4.3 Numerical simulation 

When v has a flnite value equations (30) and (31) are solved numerically. It emerges that in the 
context of the concerned scenario the universe expansion may be decelerated or accelerated from 
the start. Nevertheless once the universe starts accelerating it cannot decelerate anymore and 
eventually falls within de Sitter phase. We regard only the solutions with initial deceleration. 
Photometric observations of apparent Type la supernovae attests that the recent cosmological 
acceleration commenced at 0.3 < 2; < 0.9 [14]. To fix time scaling, a is assumed to be in the 
minimum when the redshift 2; = 0.5. 

Equations (30) and (31) are solved with initial conditions a(0) = 1 and p(0) is a constant 
greater than 1. The coefficient k is imposed to be more than •\/p(0)^ — 1 to ensure that both 
DE and DM contents are available . At the final stage DM dies out gradually. This means that 
all particles pass into the ground state and the normalized condensate number of particles in 
a comoving volume a^, namely — rica^, goes to 1 while the number density scalar ric decays 
to zero (see Fig. 1). If we start with a small value of the background number density then the 
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beginning condensate production rate is so high that not only the condensate particle number 
Nc increases but the number density ric as well. To the contrary, an excess of the condensate 
particles involves monotone behavior of the condensate number density and the particle number 
in the comoving volume. 

Fig. 2 depicts an evolution of the normalized energy densities and Q^^ of DE and DM 
respectively. The curves are plotted for different values of u and demonstrate increasing of the 
DM content when u increases. Correspondence with the current observational value of the DE 
fraction ~ 0.72 falls on u ^ 25. This implies a high-degree temperature dependence of the 
phonon energy and the second sound speed with a value much less than it would be expected 
for real phonons. Note, that in superfluid helium a lower second sound speed is provided by 
quasiparticles from the nonlinear part of the energy spectrum (such as rotons). In the context of 
the pure phonon consideration they are not taken into account and their influence is simulated 
with a large value of u. 

Plots of equation-of-state parameters wa = Pa/ Pa (A=c, n, tot) as a function of the redshift 
z is shown on Fig. 3. The DE equation-of-state Wc occurs close to —1 extremely fast after the 
universe starts to accelerate. This is accompanied by vanishing of the excitations and the total 
equation-of-state tutot is found to be the same at the same time. Simultaneously the DM one 
Wn approaches the value oi l/v that is the asymptotic value for the second sound speed as the 
sound speed Cg — > 1. If is large then DM is perceived as a nonrelativistic one. This CDM-like 
behavior is observed within the whole range of the redshift under consideration. It is expected 
that the excitations with a nonlinear dispersion will promote decreasing the equation-of-state 
Wn when the complete quasiparticle spectrum is taken into account. 

5 Conclusion 

In this letter we examine the model of superfluid Chaplygin gas (SCG) describing the dark 
sector of the universe as a matter that behaves as DE while it is in the ground state and as 
DM when it is in the excited state. Cosmological dynamics is described in the framework of 
the two-fluid model therefore the interaction between DE and DM is implicitly involved into 
the conservation laws (28) and (29). In this approach there is no need to introduce different 
equations for the description of DE and DM evolution and to use the interaction factor as an 
additional parameter. Moreover, if we abandon the self-dependent condensate ansatz (25), it 
will be impossible to uniquely divide the total energy density into the DE and DM fractions. 

The SCG model is apphed to the universe evolution from the deceleration-acceleration 
transition epoch. It provides the current mixture of the dark contents and approaches de Sitter 
phase in the future. The normal component (DM) is formed by a pure phonon gas with the 
pressure varying as the power z/ + 1 of temperature. Simple fitting the model parameters to 
the observational data shows that u must be quite large. In this case DM behaves as CDM. In 
the standard model CDM consists of nonrelativistic massive particles whereas in this model it 
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is described as an excited state of an exotic matter with a low second sound speed. To develop 
more realistic model, a wide quasiparticle spectrum should be taken into account so that the 
quasiparticles with a nonlinear dispersion are dominant in the modern epoch while the phonons 
prevail in the future at the temperatures close to zero. Further all the particles fall into the 
ground state and the final epoch is de Sitter universe. The model is integrated in more early 
epochs as a usual matter since it almost fully is in the excited state that time. 

It should be emphasized the difference between the GCG and SCG models. In the former 
the energy density ptot = Pch consists of both vacuum and matter contributions. The parameter 
a (see introduction) is assumed to satisfy the constraint a < 0.4 [17] to be in agreement with 
observations. In the latter ptot = Pc + Pn, where the condensate with the density pc = pch can be 
considered as GCG and the previous restriction is unnecessary. Moreover ordinary Chaplygin 
gas (a = 1) is preferred since a has no such a dramatic effect for this model as in [7, 8, 9] and 
its governing role for the interaction between DE and DM goes to u. This is not to say that a 
and u have an identical significance. We can see that while the GCG model becomes equivalent 
to ACDM for a = 0, SCG persists as a model unified DE and DM for i/ — > oo. Therefore we 
can use a as an additional parameter in the SCG model later. 

There is another problem in the models of DE and DM unification. It concerns the non- 
negligible sound speed that produces unphysical oscillations and an exponential blow-up in the 
DM power spectrum at present [18]. This problem was solved for the GCG model in [9] by 
the special decomposition of the energy density into DE and DM components. In the SCG 
model DM is generated by excitations and a role of the sound is played by the second sound. 
We should expect that a perturbative analysis of the energy density fiuctuations of DM will 
lead to additional restrictions on the parameters. The heuristic estimation of a relying on 
scales of Galaxy clusters obtained in [18] gives the range \a\ < 10~^ for GCG that means its 
indistinguishability from ACDM. The similar estimation for SCG can be given by the formal 
replacement a ^ 1/(1 + u). As one would expect pressureless DM {u — > oo) is free from the 
blow-up and the restriction i/ > 10^ is in agreement with the foregoing inference regarding the 
quasiparticle spectrum. 
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Figure 1: The condensate number density 
in the units of function of the red- 

shift z for k = 15 and u = 3. The sohd 
(p(0) = 1.02) and dashed (p(0) = 5) curves 
correspond to deficit and exceed of the par- 
ticles in the ground state respectively 
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Figure 2: The ratio of the energy density 
to the critical density for the different com- 
ponents as a function of the redshift z for 
k = 120 and p(0) = 1.2 
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